Effective description of hopping transport in granular metals. 
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We develop a theory of a variable range hopping transport in granular conductors based on the 
sequential electron tunnelling through many grains in the presence of the strong Coulomb interaction. 
The processes of quantum tunnelling of real electrons are represented as trajectories (world lines) of 
charged classical particles in ci+ 1 dimensions. We apply the developed technique to investigate the 
hopping conductivity of granular systems in the regime of small tunneling conductances between 
the grains g <^ 1. 

PACS numbers: 73.23Hk, 73.22Lp, 71.30.+h 
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Investigations of transport in granular conductors are 
instrumental to further advance in understanding the dis- 
ordered state of matter and became a mainstream topic 
of the current research in condensed matter physics. The 
theoretical efforts 0; IS iSi] have been triggered by 
the experiment on the granular films 0, Q, j that posed 
two fundamental questions: (i) the mechanisms of con- 
ductivity in the metallic regime where the tunneling con- 
ductance between the grains is large, g > I, and (ii) 
the origin of the exponential behavior in the insulating 
regime 
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with ctq being the high temperature conductivity that re- 
sembles Efros - Shklovskii hopping conductivity in doped 
semiconductors 0| . While the behavior of the metal- 
lic domain is now well understood P, Q , the tunneling 
transport in the insulating regime is still waiting for a 
quantitative description. The Efros - Shklovskii hopping 
conductivity in semiconductors results from the interplay 
between the exponentially falling probability of tunnel- 
ing between states close to the Fermi level, /i, and the 
Coulomb blockade effect suppressing the finite density of 
states 01 near (i.e. appearance of the soft gap Q in 
the electron spectrum). The problem of hopping trans- 
port in granular conductors is thus two-fold: (i) to ex- 
plain the origin of the finite density of states near the 
Fermi-level and (ii) uncover the mechanism of tunneling 
through the dense array of metallic grains. 

To resolve the first problem one notices that the in- 
sulating matrix in granular conductors that is typically 
formed by the amorphous oxide, may contain a deep tail 
of localized states due to carrier traps. The traps with 
energies lower than the Fermi level are charged and in- 
duce the potential of the order of e^/nr on the closest 
granule, where k is the dielectric constant of the insula- 
tor and r is the distance from the granule to the trap. 
This compares with the Coulomb blockade energies due 
to charging metallic granules during the transport pro- 
cess. This mechanism was considered in In very small 
granules one can also expect surface effects to play a role. 



Finally, in the 2D granular arrays a presence of substrate 
can also lead to the random potential. Thus in metallic 
granules the role of the finite density of impurity states 
is taken up by random shifts in the granule Fermi lev- 
els due to trap states within the insulator separating the 
granules. We describe these shifts as the external ran- 
dom potential /i;, where i is the grain index, applied to 
each site of a granular system. 

The problem of tunneling over the distances well ex- 
ceeding the average granule size [the condition necessary 
for the variable range hopping (VRH) mechanism to oc- 
cur] in a dense granular array has not been resolved: no 
mechanism by which the sequential tunnelling in a sys- 
tem of many grains can proceed were known neither the 
technique that could provide the proper description of 
such processes in the presence of strong Coulomb inter- 
actions were available. The virtual electron tunnelling 
through a single granule (quantum dot) was considered 
in Ref. T?]. However, the method of [3| does not allow 
straightforward generalization on the granular arrays. A 
phase functional approach introduced in Ref. that was 
recently applied to the description of metallic granular 
samples at not very low temperature, T > gS, [ij, where 
S is the average mean energy level spacing for a single 
grain, could have been considered as a possible candidate 
for such a description. However, this approach includes 
only nearest-neighbor electron tunneling processes and 
thus is not capable of capturing the features of coherent 
sequential electron tunneling through many grains. 

In this Letter we develop a technique providing a sim- 
ple and transparent description of the sequential virtual 
tunnelling through the granular arrays and offering a 
ground for a comprehensive theory of the hopping trans- 
port in granular metals. Our method is based on two 
ideas: (i) following the approach of Ref. we absorb 
Coulomb interactions into the phase field and (ii) using 
the method of Ref. we map the quantum problem that 
involves functional integration over the phase fields onto 
the classical model that has an extra time dimension. 
The advantage of such a procedure is that the higher 
order tunnelling processes can be included in a straight- 



forward way along with the nearest neighbor hopping. 

Using the developed technique we find the probabil- 
ity of a sequential tunnelling through several grains. For 
the diagonal Coulomb (short range) interaction the tun- 
nelling probability between the 0th and Kth grain can 
be written as a product of the sequential probabilities Pk 
as 
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where gk is the tunneling conductance between the fcth 
and [k + l)st grain and E"^ = E^ ± with E^ being 
the charging energy of the fcth grain. Upon averaging, the 
hopping probability falls exponentially with the distance, 
P e~^/^, where the localization length ^ is 



f ^ ca/\ii{ETr/g6). 



(3a) 



Here a is the grain size, c = 1 is the numerical constant, 
E and g represent geometrical averages of the energies, 
Ek = 2 {^/E^ + 1/E^)~^, and tunnelHng conductances 
respectively taken along the optimal path corresponding 
the maximal probability of tunneling between the initial 
and final granules: 
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Note that the optimal path according to Eqs. © is de- 
termined by contributions from both conductances and 
effective Coulomb energies Ek- Result (I3a|l holds also for 
the long range Coulomb interactions which only renor- 
malize the numerical constant to some 0.5 < c < 1. 

Taking a typical value for the tunneling conductance 
gk — g — iV exp(— c?/A), where N is the number of chan- 
nels in a tunnelling contact, d is the thickness of the 
insulating layer and A is the localization length within 
the insulator, the localization length ^ in Eq. I|3a(l can be 
estimated as - ca [d/X + ln{E/NS)] ~\ 0- 

Applying conventional Mott-Efros-Shklovskii argu- 
ments and using Eq. (|3a|l we obtain Eq. ^ for the hop- 
ping conductivity with the characteristic temperature 
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K is the effective dielectric constant of a granular sample. 

Now we outline our approach and derive the formula 
(|3a|l . which is one of our main results. We consider a d- 
dimensional array of metallic grains. Electrons move dif- 
fusively inside each grain and tunnel from grain to grain. 
We assume that the dimensionless tunneling conductance 
g is smaller than the intra-granule conductance. The sys- 
tem of weakly coupled metallic grains is described by the 
Hamiltonian 



H — Hq + He 
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FIG. 1: Diagram a) represents the lowest order correction 
in tunneling conductance g to the partition function, Eq. @ . 
The solid lines denote the propagator of electrons, the crossed 
circles stand for the phase dependent tunnelling matrix ele- 
ments tije*'^'^''^'. Diagram b) describes a higher order correc- 
tion in tunnelling matrix elements tij that includes impurity 
averaging within each grain. The dotted lines represent im- 
purity scattering. The diagram c) shows the zero dimensional 
diffuson in time representation. 



where ty is the tunneling matrix element corresponding 
to the points of contact ri and Vj of ith and jth grains. 
The Hamiltonian Hq in the r. h. s. of Eq. (|5a|) describes 
noninteracting isolated disordered grains. The term He 
includes the electron Coulomb interaction and the local 
external electrostatic potential /ii on each grain 



He = 



(5b) 



where Cy is the capacitance matrix and fii is the oper- 
ator of electron number in the ith grain. The Coulomb 
interaction term written in a Lagrangean form can be 
decoupled with the help of the potential field Viir): 

^c = --^'^[Vi+ i^ii] Cij [Vj + -i'^riiVi, (6) 



where rii ~ ni{T) is the electron density field. The last 
term in the r. h. s. of Eq. ^ can be absorbed by the 
fermion gauge transformation ^iij) — > ipiiT) e~ 



I Mr) 



Mr) 



V,{T')dT', 



such that charge action becomes 

^c = -^YJ^V^+ i^lr] C^j [Vj + iflj] 



(7) 



(8) 



Since the action of an isolated grain is gauge invariant the 
phases 4>i (r) enter the full Lagrangean of the system only 



through the tunnelling matrix elements ti 



where (jjijir) — (j)i{T) — (f>j{T) is the phase difference of the 
ith and jth grains. We will construct the effective charg- 
ing functional using the diagrammatic expansion of the 
partition function Z in tunnelling matrix elements |l6l |. 
The lowest order correction Zi to the partition function 
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FIG. 2: Illustration of the mapping of a quantum granular 
array model onto the classical model of the charges imbedded 
into the electron world lines forming loops. Figure (a) de- 
scribes a quadruple while figure (b) represents the K-th order 
loop. Though the illustration is shown for the case of 1 + 1 
dimensions (Id granular array) generalization to higher d+1 
dimensions is straightforward. 

Z defined as Z = Zo{l + Zi + ...) with Zq being the 
partition function of the isolated grains is given by the 
diagram shown in Fig. la. The correction Zi at zero 
temperature is given by the expression 

^1 ['^] = ^ L 5y J dn dT2 _ , (9) 

yet to be averaged over the phases (f>i{T). The discon- 
nected diagrams must also be included in the partition 
function such that the complete partition function that 
takes into account nearest neighbor electron tunneling 
becomes 

oc 

^jvjv = ^o 51(^1 [<^])/fc!- (10) 
fc=i 

Here the angular brackets stand for the averaging over 
the phase fields (j)i{T) with the charging action ((SJ. 
Equation (|IO|l coincides with the partition function ob- 
tained from the phase functional of Refs.fllll^bv expan- 
sion in small tunneling conductance g. Integration over 
the phase fields 4>i{T) can be implemented exactly for 
each term in the expansion H10|l using Eqs. © and 
with the help of the formula 

^e«i:„0.(T„)e„\^ ^g-t/^ (11) 

where U can be viewed as a Coulomb energy of a sys- 
tem of classical charges in a space with one extra time 
dimension interacting via the Id Coulomb potential 0| 

f/ = ^ ^ -Ec"'*"^|T„j-r„Je„i e„2+^A^i„T„e„. (12) 

ni,n2 n 

The physical interpretation of Eq. (|12|) is as follows: the 
charges e„ = ±1 interact via the Id Coulomb potential 



along the time direction while the strength of the inter- 
action is given by the charging energy Ec = e^C~^ jl. 
In addition, there is a linear growing potential \XiT that 
can be viewed as a local constant electric field applied 
in a time direction. 

The partition function, Eq. (|10l) , of the quantum prob- 
lem under consideration in the nearest neighbor electron 
tunneling approximation can be presented as a partition 
function of classical charges that appear in quadruples 
as shown in Fig. 2a. Each quadruple has its "internal" 
energy fn^^ where the upper index stands for quadru- 
ple approximation, that comes from the electron Green 
function lines in Fig. la 

4i)=ln[2^(AT„)Vg„]. (13) 

Here At„ is the size of the nth quadruple in the r— direc- 
tion in Fig. 2a, gn is the tunneling conductance between 
the sites occupied by the quadruple and all charges are 
subject to the Coulomb interaction and local potentials 
in accordance with Eq. ()12|l . Thus the total energy of a 
system of N quadruples can be written as 

N 

E^'^=J2£k'^+U4N, (14) 

n=l 

where Sn"^ is the internal energy of the nth quadruple and 
UiN is the static Coulomb energy of 4N charges that form 
N quadruples defined by Eq. H12() . We notice that the dif- 
ferent quadruples interact with each other only through 
the Coulomb part of the energy U. The model (|14f) can be 
understood as a classical analog of the phase functional 
of Rcf. 13 and thus it has the same region of applicability. 
Wc see that within this approximation the hopping con- 
ductivity obviously cannot be described and that higher 
order tunnelling events must be included. 

The advantage of mapping the original quantum model 
on the classical electrostatic system is that it allows 
to include the higher order sequential tunneling pro- 
cesses shown in Fig. 2b in essentially the same way as 
the nearest neighbor hopping. The higher order dia- 
grams include the single grain diffusion propagator Dq — 
2TrS~^ /\oJi — UJ2\ that having being transformed into the 
time representation result in the expression shown in 
Fig. Ic. Each tunneling matrix element contains phase 
variables as tije^'^^^'^'^\ Thus, the diagram in Fig. lb can 
be presented as a charge loop shown in Fig. 2b. The 
charges interact electrostatically and are subject to the 
local potentials according to Eq. (|12|l . The internal en- 
ergy of a single K—th order loop is 

= ln[27ri^/go] - XI +l^l^^ol 
fc=i 

K~l 

+ ln|Arx|+Xln(|Arfe| + |Ar^|), (15) 

k=l 
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where the time intervals Ar and Ar' shown in Fig. 2b 
have opposite signs. 

Now we show how the effective classical model can be 
used to derive the tunnelling probability through several 
grains that is the key quantity for understanding the hop- 
ping conductivity regime. Let us consider two sites where 
the local Coulomb gap is absent. In our model it means 
that the Coulomb energy on these two sites is compen- 
sated by the external local potential ^i. Thus, removing 
the charge from the grain iq as well as placing the charge 
on the grain does not cost any Coulomb energy. The 
electron loops that appear in the classical representation 
can be viewed as electron worldliness. Thus the proba- 
bility of tunneling between two sites is given by the free 
energy of the configuration shown in Fig. 2b. The lower 
part in Fig. 2b represents the tunneling amplitude A from 
the site io to the site while the upper part represents 
the inverse process. The probability P — A* A is given by 
the whole loop while the internal energy corresponding 
to the initial and final states should not be counted. 

The contribution of the loop of the K-th order shown in 
Fig. 2b can be calculated easily for the case of a diagonal 
Coulomb potential = SijE^,. In this case integrations 
over the time intervals Ar/j and Ar^ can be done inde- 
pendently and the tunneling probability P is given by 
the product of sequential probabilities 

. M / (16) 

resulting in expression for Pk given by Eq. jSJ. 

In the presence of the long range Coulomb interactions 
the situation is more complicated since the integrals over 
variables Ar^ , At^ cannot be taken on each site k in- 
dependently. However one can estimate the tunnelling 
probability by finding its upper and lower limits: Indeed, 
certainly just neglecting the long range part of the po- 
tential (as we did above) one gets the upper boundary. 
The lower boundary can be obtained by considering such 
trajectories where all At^ are of the same sign. In such a 
case the electric field created by the charges on step k do 
not interfere with other charges and thus one can imple- 
ment integrations independently. This will give a lower 
boundary for the tunnelling probability which is about 
a factor 2 smaller than the tunnelling probability Pk ob- 
tained neglecting the long range part of the Coulomb po- 
tential. This results in the boundaries of the pre factor 
0.5 < c < 1.0. 

The above derivation of tunnelling amplitudes was 
done for the zero temperatures case (1^. Of course the 
finite temperatures and the presence of phonons are nec- 
essary for the realization of the Mott-Efros-Shklovskii 
mechanism. While there is a temperature interval where 
temperature effects appear only as a pre-exponential fac- 
tor and do not interfere to tunnelling probabilities, the 
extension of the developed method to finite tempera- 
tures remains an important task. At zero temperature 



the system occupies the unique ground state, whereas at 
finite temperatures the exited states must be included as 
well. Such exited states can be described by the introduc- 
tion of winding numbers representing (after Poisson re- 
summation) the static charges in the system. This gen- 
eralization however goes beyond the scope of the present 
Letter and will be a subject of forthcoming publication. 

In conclusion, we have developed the technique en- 
abling quantitative description of hopping conductivity 
of granular conductors in a low, g ^ 1, tunneling conduc- 
tance regime and derived the Coulomb blockade-governed 
VRH conductivity a oc exp[— (To/T)^/^] of granular ma- 
terials. The essential feature of our approach is represen- 
tation of sequential intergranular quantum tunneling of 
the electrons as trajectories of charged classical particles 
in a + 1 dimensional system. 
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